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The Binder cumulant at the phase transition of Ising models on square lattices with ferromagnetic 
couplings between nearest neighbors and with competing antiferromagnetic couplings between next- 
nearest neighbors, along only one diagonal, is determined using Monte Carlo techniques. In the phase 
diagram a disorder line occurs separating regions with monotonically decaying and with oscillatory 
spin-spin correlations. Findings on the variation of the critical cumulant with the ratio of the 
two interaction strengths are compared to related recent results based on renormalization group 
calculations. 
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In the field of phase transitions and critical phenom- 
ena, the fourth order cumulant of the order parameter—, 
the Binder cumulant U, plays an important role. In par- 
ticular, the cumulant may be used to locate the phase 
transition from the intersection of the cumulant for dif- 
ferent system sizes. The cumulant also allows to compute 
the critical exponent of the correlation length, and thence 
to identify the universality class of the transition. 

However, care is needed when attempting to identify 
the universality class from the value of the Binder cumu- 
lant, in the thermodynamic limit, at the transition, U*. 
Indeed, that value is known to depend, in a given univer- 
sality class, on various aspects including the boundary 
conditions, the shape of the system (being extrapolated 
to the thermodynamic limit), and the anisotropy of the 
correlations or interactionsjii^i^i^i^i^ii On the other hand, 
U* may not depend on other details of the system like 
the spin valued or the lattice structure 

In recent years, based on renormalization group cal- 
culations of Dohm and Chen^i^ and subsequent Monte 
Carlo simulations^i^iiS, the influence of anisotropic in- 
teractions on the critical Binder cumulant U* has been 
elucidated. Much attention has been focused on the Ising 
model with nearest neighbor (nn) couplings, J, where the 
anisotropy is introduced by the next-nearest (nnn) cou- 
plings along only one diagonal of the lattice, Jdr^- One 
encounters a 'nondiagonal anisotropy matrix—. Then, in 
general, there is no simple transformation relating [/*, 
for given boundary condition and shape, to that of the 
isotropic model by adjusting the shaped. Such a tran- 
scription may be easily performed in the case of a diag- 
onal anisotropy matrix, as it occurs, for instance, for the 
nn Ising on a square lattice with different vertical and 
horizontal ferromagnetic interactions, where U* of the 
anisotropic model on lattices with square shape may be 
expressed by U* of the isotropic model on lattices with 
rectangular shapes.— 

In this contribution, we shall extend our previous 
Monte Carlo study for anisotropic nn and nnn Ising mod- 



els on square lattices with only ferromagnetic interac- 
tions, J, Jd > 0, to the case of competing nnn antiferro- 
magnetic couplings, Jd < 0, with J remaining ferromag- 
netic. Thereby, spatially modulated, oscillatory spin- 
spin correlations may occur, adding an interesting fea- 
ture to the phase diagram. The present study has been 
partly motivated by related recent quantitative renormal- 
ization group calculations of U*{Jd/J) for a closely re- 
lated model, showing intriguing symmetry properties^, 
as will be discussed below. 

The Hamiltonian of the model may be written in the 
form 



H — — Sx,y{J {Sx+l.y + Sx,y+l) + JdSx + l,y + l)) (1) 

x,y 

where Sx.y = ±1 is the Ising spin at site {x,y). J > is 
the ferromagnetic nn coupling along the principal axes of 
the square lattice, while the nnn coupling, Jd acts along 
only one diagonal, i.e. along the [11] direction of the 
lattice. Jd may be ferromagnetic, Jd > 0, as has been 
studied before^, or antiferromagnetic, Jd < 0. 

Before turning to the critical Binder cumulant, we shall 
first discuss the phase diagram, exhibiting interesting fea- 
tures, especially for antiferromagnetic nnn couplings. 

At Jd/J > —1, there is a ferromagnetic phase at 
low temperatures. The exact transition temperatures, 
ksTc/J, are known to be determined b y^^'^^ 

{smh(2 J / keTc) f +2 smh{2 J/ kBT,.)smh{2Jd/kBTa) = 1, 

(2) 

The line is depicted in Fig. 1. The transition temper- 
ature Tc goes to zero on approach to J^/J = —1. At 
that point, the ground state is highly degenerate, with 
the energy per site being Eq— — J. The degenerate con- 
figurations include the ferromagnetic structures, uncou- 
pled antiferromagnetic Ising chains along the [11] direc- 
tion, and horizontal and vertical ferromagnetic stripes 
of spins with alternating sign corresponding to coupled 
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modulated Ising chains along the [11] direction. A rather 
large degeneracy, due to the uncoupled antiferromagnetic 
chains along the [11] direction, persists at Jd/J < — 1- 
Indeed, it has been argued that there is no long-range 
ordering at low temperatures in that part of the phase 
diagranJ^. Actually, in the following we shall consider 
Jd/J > -1. 
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to the [11] direction, G3(r)= < Sx,ySx+r,y-r >■ At 
T > Td, Gi and G2 decay, for sufficiently large distances 
r, exponentially and in an oscillatory, purely sinusoidal 
manner, with the wavenumber depending on Jd/J and 
temperature^^. At and below the disorder line, the cor- 
relations along the principal axes and along the [11] di- 
rection decay monotonically. In addition, our simulations 
suggest that G3 decays also above Td monotonically (and 
exponentially) with distance, indicating that the compet- 
ing interactions do not affect qualitatively correlations 
perpendicular to the [11] direction. 
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FIG. 1: Phase diagram of Hamiltonian (1), showing the exact 
boundary line of the ferromagnetic phase^^ (solid line) and 
the disorder line^^ (dashed line). 

Due to the competing ferro- and antiferromagnetic 
couplings, at Jd/J < 0, a disorder hue, Td{Jd/J), sepa- 
rates, above Tc, the region with only monotonically de- 
caying spin-spin correlations from the one with oscilla- 
tory correlational^. The disorder line of the model, eq. 
(1), has been calculated exactly as wel h^^'^^ It is deter- 
mined by 

cosh(2 J/fcsTd) = exp(-2Jrf/fcsTd). (3) 

The line Td{Jd/J) is also shown in Fig. 1. It arises from 
the highly degenerate point at {Jd/J = —1,7" = 0), re- 
flecting the spatially modulated configurations occurring 
there as ground states. At low temperatures, ~Jd be- 
ing not far from J, the disorder line Td follows closely 
the phase transition line, Tc, eventually moving away 
from it towards higher temperatures, as Jd gets weaker, 
and finally approaching infinite temperature at vanishing 
nnn antiferromagnetic couplings, see Fig. 1. Of course, 
there is no disorder line for ferromagnetic nnn interac- 
tions, Jd > 0- 

The different types of correlations at temperatures be- 
low and above the disorder line, fixing Jd/J, are il- 
lustrated in Fig. 2, displaying Monte Carlo data for 
spin-spin correlations along the [11] direction, Gi(r)= 

< Sx.ySx+r.y+r >, along the principal axes, G2(r)= 

< Sx,ySx,y+r >=< Sx,ySx+r,y >, and perpendicular 



FIG. 2: Simulated spin-spin correlation function Gi, with i= 
1 (circles), 2 (squares), 3 triangles), as described in the text, 
at Jd/J ~ —3/4 and ksT/ J = 0.9, for lattices with 48^ spins. 

Note that disorder lines also exist in other Ising models 
with competing interactions, for instance, in the much 
studied ANNNI modelJ^Ji 

Let us now turn to the Monte Carlo findings on the 
critical Binder cumulant U* of the model, eq. (1). U* is 
defined byi 

U* = C/(Tc) = 1- < > /(3 < Af2 >2) (4) 

taking the thermodynamic limit; < > and < M* > 
denote the second and fourth moments of the order pa- 
rameter, the magnetization M. 

To estimate U*{Jd/J), we simulated the model for 
square shapes with sites or spins, employing full peri- 
odic boundary conditions, using the standard Metropo- 
lis algorithm (note that, e.g., cluster ffip algorithms are 
usually rather inefficient in case of competing ferro- and 
antiferromagnetic interactions). Monte Carlo runs with, 
typically, 5 x 10^ Monte Carlo steps per site were per- 
formed, averaging then over several, up to ten, of these 
runs to obtain final estimates and to determine statisti- 
cal error bars. L ranged from 4 to 64. To extrapolate 
to the thermodynamic limit, L — > 00, least square fits 
were done. The procedure is exemplified in Fig. 3. The 
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FIG. 3: Binder cumulant at Tc, U{Tc), versus inverse sys- 
tem size 1/L, including least square fits to the Monte Carlo 
data (dashed lines) and the resulting estimates for the criti- 
cal Binder cumulant, U* . The cases Jd/J= —0.316 (squares) 
and —5/9 (circles) are shown. 

final error bars result from the fits (one may emphasize 
that the finite-size behavior of the Binder cumulant is 
not known in the anisotropic case). 

In Fig. 4, we plot U* against s = 1/(1 + {J/Jd)), 
following a recent renormalization group analysis^., with 
s < for < and s > for > 0. The reason 
for this choice will be discussed below. We also include 
our previous estimates^, of U* for Jd > 0, adding one 
additional new point at s = 1/6. Obviously, the criti- 
cal Binder cumulant depends, both in the case of only 
ferromagnetic interactions and also in the case of com- 
peting antiferromagnetic couplings, continuously on s (or 
Jd/ J)- The transition at Tc{s) is always in the Ising uni- 
versality class, as we confirmed by monitoring, especially, 
the, asymptotically, logarithmic size dependence of the 
specific heat at Tc. Thence, U* , for given boundary con- 
dition and shape, changes, within the Ising universality 
class, with the anisotropy, Jd/ J- 

Interestingly, U*{s) is (almost) symmetric around the 
isotropic situation, s = 0, at small values of |s|. However, 
for larger values of |s|, there is a pronounced asymme- 
try. This may be understood by the observation that the 
model reaches one-dimensional limits, where U* — > 0, 
for s = 1 or Jd/J = oo at positive s, and for s — —oo 
or Jd/ J = — 1 at negative s. The nonmonotonicity of U* 
for fairly small values of \s\ is also worth mentioning. It 
has been first found for positive s in a previous Monte 
Carlo studyS., and it also holds for negative s. 

The previous Monte Carlo data^ on U*{s) at s > 
agree closely with results of a recent renormalization 
group (RNG) analysis by Dohm, as has been already 
discussed there-. That analysis has been done for the 
model, eq. (1), augmented by a ferromagnetic interlayer 
nn coupling Ji on a cubic lattice, with Ji = J + Jd- 



FIG. 4: Critical Binder cumulant, U*{s), versus s = 1/(1 + 
{J / Jd))- The result of a renormalization group analysis^, see 
text, is included as dashed line. 



The interlayer coupling has been chosen to preserve es- 
sential features of the anisotropy matrix of the model on 
the square latticeji Indeed, this choice allows to describe 
well, in the framework of RNG theory, the simulational 
data on U*{s) — U*{0) for the two-dimensional model 
with ferromagnetic nnn interactions, s > 0.~ Certainly, 
the critical Binder cumulant in the isotropic case, s = 0, 
is different for square and cubic lattices, thereby motivat- 
ing to look at the difference, J7*(s) — J7*(0), of the critical 
cumulants^. The RNG study yields a perfect symmetry 
of U*{s) — J7*(0) around s = 0, and an interesting non- 
monotonic behavior in s. 

To compare the present Monte Carlo data to the RNG 
findings, we prefer to plot, see Fig. 4, the RNG re- 
sults in the form U^jyQ{s)U*{s = 0)/[/^^q(s = 0), 
with the critical Binder cumulant of the two-dimensional 
isotropic Ising model, U*{0), being^ 0.61069... Of course, 
the perfect symmetry and nonmonotonicity found in the 
RNG study is preserved by this choice, as displayed in 
Fig. 4. As follows from that figure, the symmetry, sug- 
gested by the RNG analysis, is approximated closely 
in the two-dimensional Ising model with antiferromag- 
netic nnn interactions for small values of |s|, confirming 
the nonmonotonicity also for negative s. But there are 
pronounced deviations at stronger antiferromagnetic nnn 
couplings. 

In fact, for competing antiferromagnetic couplings the 
phase diagrams of the models in two and three dimen- 
sions, and thence the corresponding critical Binder cu- 
mulants, may be expected to differ substantially. In par- 
ticular, in two dimensions there is, as discussed above, 
a phase transition of Ising type down to Jd/J — — 1, 
with spatially modulated correlations occurring only in 
the disordered phase above the disorder line. In con- 
trast, in three dimensions, the transition line of Ising 
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type between the ferro- and paramagnetic phases may 
be expected to extend only up to a Lifshitz poin t"'^^'^^i^° 
at {Jd/ J)lp (> —1) or Sip, at which the ferromagnetic, 
the spatially long-range ordered modulated, and the dis- 
ordered phases meet, similar to the well-known situation 
in the ANNNI model^ 

Indeed, a standard mean-field calculation shows that 
the Lifshitz point occurs, for Ji = J + Jd, at Jd/J = 
— 1/2, corresponding to s^p = —1. In the modu- 
lated phase close to the transition to the paramagnetic 
phase, at Jd/J < {Jd/J)LP = ^1/2, the magnetization 
along the x- and y- axes change in a sinusoidal man- 
ner. The modulated, long-range ordered phase seems 
to arise from the highly degenerate ground states at 
Jd/J = —1, with the transition line between the ferro- 
magnetic and the modulated phases, at — 1 < s < —1/2, 
being of first order. The vicinity of that degenerate point 
may be explored by systematic low-temperature series 
expansions'^, being however, well beyond the scope of the 
present study. Because there are modulated magnetiza- 
tion pattern along two principal axes of the cubic lattice, 
the X- and the j/-axes, we are dealing here with a 'biaxial 
Lifshitz point' (m 2 in the standard notatior^i^). 

The existence of a Lifshitz point at s = —1, has been 
argued to allow for the complete symmetry of U*{s) 
around s = in the RNG analysis^. Going beyond mean- 
field theory, the location of the biaxial Lifshitz point may 
shift somewhat. To determine accurately the position of 
the Lifshitz point, high temperature series may be very 
helpful, as has been found, for instance, in the case of 
the ANNNI and related models i^i^^i^^ Obviously, it is 
an open problem, in which way the possible shift in the 
Lifshitz point may affect the proposed perfect symmetry 
of U* (s) . Note that the lower critical dimension of a bi- 
axial Lifshitz point is thre o^^'^"^ , excluding its existence 
at non-zero temperatures for the square lattice, but not 



for the cubic lattice. 

Certainly, it is desirable to determine U*{s) of the 
three-dimensional model in simulations. However, nei- 
ther Tc{s) nor the location of the Lifshitz point are 
known. To get then data on U* of the required high accu- 
racy, a huge amount of computer time would be needed, 
being hardly feasible at present. 

In summary, the critical Binder cumulant U* of an 
anisotropic two-dimensional Ising model with compet- 
ing ferro- and antiferromagnetic interactions has been 
determined. Because the transition is known exactly, 
one arrives at accurate estimates based on extensive 
Monte Carlo simulations. Employing full periodic 
boundary conditions and considering square shapes, U* 
is found to vary continuously with changing anisotropy 
s. A remarkably close agreement of our findings on 
[/*(s), at positive values of s and, in case of competing 
interactiosn, for small negative values of s with the 
results of a recent renormalization group analysis is 
observed. Differences at stronger competing anisotropy 
are explained by the absence of a biaxial Lifshitz point 
at non-zero temperatures in two dimensions. 
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